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Radar Precoder Design for Spectral Coexistence 
with Coordinated Multi-point (CoMP) System 

Jasmin A. Mahal, Awais Khawar, Ahmed Abdelhadi, and T. Charles Clancy 


Abstract —This paper details the design of precoders for 
a MIMO radar spectrally coexistent with a MIMO cellular 
network. We focus on a coordinated multi-point (CoMP) sys¬ 
tem where a cluster of base stations (BSs) coordinate their 
transmissions to the intended user. The radar operates in two 
modes, interference-mitigation mode when it avoids interference 
with the CoMP system and cooperation mode when it exchanges 
information with it. Using either the conventional Switched Null 
Space Projection (SNSP) or the newly proposed Switched Small 
Singular Value Space Projection (SSSVSP), the radar beam 
sweeps across the BS clusters focusing on the optimal ones, 
optimal in either nullity or difference between the precoded and 
original radar signal. Taking the channel estimation error into 
account, the design of precoder is pivoted on the minimal radar 
interference at the BS clusters during Interference-mitigation 
mode and minimal bit-error-rate at the BSs during cooperation 
mode while interfering minimally with the radar target detection 
capability. Our investigation shows that loss in radar performance 
can be compensated using SSSVSP instead of SNSP to some 
extent but increasing the number of radar antenna elements 
goes a long way to serve the purpose. Simulations verify our 
theoretical predictions about the proposed SSSVSP. 

Index Terms —MIMO Radar, Null Space Projection, Spectral 
Coexistence, CoMP, Spectrum Sharing, MMSE, Zero Forcing, 
Linear Precoding 


I. Introduction 

Cellular network operators are predicting a 1000 x increase 
in capacity to keep pace with the tremendous growth of 
aggregate and per-subscriber data traffic III. Increased 
investments in infrastructure, e.g. larger number of small 
cells, and more spectrally efficient technologies, e.g. LTE- 
Advanced, can help meet this challenge partially. The 
Federal Communications Commission (FCC) is considering 
a number of options including incentive auctioning and 
sharing of federal spectrum to meet the commercial spectrum 
requirements. Of these two, spectrum sharing is a quite 
promising technology due to the large number of under¬ 
utilized federal spectrum bands that can be shared with 
commercial cellular operators to satisfy their growing 
demands. But spectrum sharing is associated with its inherent 
set of challenges because the incumbents need to be protected 
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from the harmful interference that can arise due to the 
operation of other systems in the shared bands. In this paper, 
the authors have addressed the specihc issue of spectrum 
sharing between a multiple-input multiple-output (MIMO) 
radar system and a commercial MIMO cellular network 
consisting of clusters of cooperative base stations (BS), 
commonly known as coordinated multi-point system. 

Coordinated Multi-point (CoMP) Systems: MIMO is 
a key enabling technique for improving the throughput of 
future wireless broadband system. Despite its theoretical 
attractiveness, a MIMO commercial cellular system is 
essentially interference-limited because the transmission in 
each cell causes interference to the other cells, i.e. inter-cell 
interference (ICI). For cell-edge user equipments (UFs), ICI 
is severe due to their proximity to a number of neighboring 
cells. Downlink is more interference-limited than the uplink 
because the sophisticated interference suppression techniques 
can be easily implemented in BSs but they are not practical 
for UFs. Implementation of these techniques would render the 
UFs power-inefficient or incompact. Moreover, coordination 
among UFs to mitigate interference is challenging and, thus, 
does not alleviate the ICI problem. Conventional approaches 
for mitigating the ICI, such as static frequency reuse, spread 
spectrum, and sectoring etc. are not efficient for MIMO 
networks owing to their respective inherent limitations ||2l. 
Due to the recent advancements in processing capabilities 
at BSs and the increase in backhaul capacity, coordination 
among BSs can be utilized to mitigate the ICI. Such a 
configuration is known as coordinated multi-point (CoMP) 
transmission/reception which not only controls the ICI but 
also takes advantage of it. Due to the advantages of the CoMP 
system over traditional cellular system it has been included 
in the fourth generation (4G) mobile standard, i.e., 3GPP 
LTE-Advanced Release 11 and beyond consider CoMP as 
an enabling technology for 4G mobile systems IH. In short, 
the CoMP system coordinates simultaneous transmissions 
from multiple BSs to UFs in the downlink and perform joint 
decoding of UE signals at multiple BSs in the uplink. This 
results in improved coverage, throughput, and efficiency for 
cellular system. 

Regulatory Efforts to Share Radar Spectrum: In 2010, 
the U.S. President Obama issued a Presidential Memorandum 
which recommended the Federal Government to make avail¬ 
able 500 MHz of Federal or non-federal spectrum for both 
mobile and fixed wireless broadband use by commercial users 
within 10 years ||4|. In compliance with the Presidential direc- 
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tive, the National Telecommunications and Information Ad¬ 
ministration (NTIA) identified a number of bands compatible 
for commercial utilization including the 3500-3650 MHz band 
which is primarily used by federal radar systems El. However, 
the NTIA also pointed out in a study that if the identified 
radar bands are shared with commercial wireless services, 
exclusion zones are required to protect wireless services and 
these exclusion zones may go upto hundreds of kilometers m. 
On the other hand, the Federal Communications Commission 
(FCC) has proposed to use small cells that will be sharing 
the 3.5 GHz band with the federal radar systems I?). In the 
past, government bands have been successfully shared with 
commercial wireless systems. Example includes; Wi-Fi and 
Bluetooth in the 2450-2490 MHz band and wireless local area 
network (WLAN) in the 5.25-5.35 and 5.47-5.725 GHz radar 
bands IH. 

A. Related Work 

Various schemes have been proposed to share radar spec¬ 
trum with communication systems. Spectrum sharing schemes 
based on waveform shaping was first proposed by Sodagari et. 
al. M in which the radar waveform was shaped in a way so 
that it was in the null space of interference channel between 
the radar and a communication system. This scheme was 
extended to multi-cell communication systems by Khawar et. 
al. ITOl in which spectrum was shared between a MIMO radar 
and a LTE cellular system with multiple base stations. They 
proposed algorithms to select the interference channel with the 
maximum null space to project the radar signal onto that null 
space. Although this paper addressed a more practical scenario 
with multiple BSs, at any given time it mitigated interference 
to only one BS, i.e, the station with the maximum nullity. 
A similar approach was presented by Babaei et. al. in which 
radar waveform was shaped to mitigate interference to all the 
BSs in the network ED. Other spectrum sharing approaches 
include cooperative sensing based spectrum sharing where a 
radar’s allocated bandwidth is shared with communication sys¬ 
tems lll2l- lfT6l . joint communication-radar platforms that are 
spectrally-agile ED, database-aided sensing at communication 
systems to enable radar spectrum sharing ED, radar systems 
that form virtual arrays to coexist with communications system 
ED, and beamforming approaches at MIMO radars ll20l . 

B. Our Contributions 

Building upon the works in lITOl and ED, the authors of this 
paper have extended the solution approach to address a MIMO 
radar and a commercial CoMP communication system coex¬ 
istence scenario, which is applicable for the LTE-Advanced 
system. Our contributions are summarized as follows: 

• Precoder Design for Interference Mitigation: In order 
to mitigate interference to CoMP system we present two 
novel radar precoder design methods in this paper. In the 
first scheme, the radar projects its signals onto the clusters 
formed by the CoMP system and selects the cluster 
with the maximum null space. This scheme is known as 
Switched Null Space Projection (SNSP) as radar looks for 
the optimal cluster at each pulse repetition interval (PRI) 


and thus switches among clusters at each PRI depending 
upon its mission requirements. In the second scheme, 
called Switched Small Singular Value Space Projection 
(SSSVSP), the projection space has been expanded to 
include the subspace corresponding to the small non-zero 
singular values under a specified threshold in addition to 
the null space. The precoder is designed based on the 
knowledge of a composite interference channel matrix 
which comprises of the channel matrices from radar to 
all the BSs in a particular cluster. Based on the proposed 
channel estimation technique in ||2D . all of the BSs in the 
cluster coordinate in their choices of training symbols and 
power transmission during the channel estimation phase 
and the radar transmitter can estimate the composite 
interference channel between them. 

• Precoder Design for Cooperation: We design a novel 
radar precoder for communicating with the CoMP system. 
The purpose of this precoder is information exchange 
between radar and the CoMP system. The CoMP system 
primarily informs the radar of the clustering information 
and in return radar informs the CoMP system about which 
cluster it has selected for interference mitigation. Basi¬ 
cally, this precoder is vital for the success of interference 
mitigation techniques. Without this sort of information 
exchange, spectrum sharing would not be successful. We 
design simple linear precoders based on zero forcing (ZE) 
and minimum mean-square-error (MMSE) criteria. 

> Precoder Performance Analysis: In order to evaluate the 
performance of radar precoders, we perform simulations 
in the presence of channel estimation errors ||22l. We 
look at the target localization and interference mitigation 
capabilities of the radar precoders. The results indicate 
that while the precoder nullifies the radar interference 
to the clusters, it degrades the radar performance by 
introducing correlation in the probing signals. We show 
that this performance loss can be compensated for by two 
means either by increasing the number of radar antennas 
or by utilizing small singular value space projection rather 
than using null space projection. Our results show that 
between the two compensation schemes, the former is 
more effective. 

C. Notations 

The vectors and matrices are denoted by lower-case and 
uppercase boldface letters, respectively (e.g., j and J). The 
rank, null space, transpose, and Hermitian transpose of J 
are represented by rankjj}, N {J}, and J* , respectively. 
The subspace spanned by a set of vectors S is denoted 
by Span {S'}. The by TV identity matrix is presented by 
I^v. 11-112 IMIf denote L-2 norm and Erobenius norm 
respectively. Eor a quick reference, important notations are 
summarized in Table |I] 

D. Organization 

The remainder of the paper is structured as follows. Section 
inidetails the radar/communication system spectral-coexistence 
model. Design of radar precoder for different operating modes 
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TABLE I 

Table of Notations 


Notation 

Description 

M 

Total number of base stations (BSs) 

K 

Total number of user equipments (UEs) 

Mfc 

Cluster of BSs serving fc* UE 

Mk 

Total number of BSs in cluster 

Mr 

Radar transmit/receive antennas 

A'bs 

BS transmit/receive antennas 

Mje 

UE transmit/receive antennas 

■frk 

^R,i 

Composite interference channel between radar 

and the cluster 


is described in Section Hin The newly proposed SSSVSP 
algorithms are presented in Section |IVl Section |V] explains 
the simulation results. Section |Vl] concludes the paper with 
final remarks. 

II. Radar/CoMP System Spectral-Coexistence 
Models 

In this section, we introduce CoMP communication sys¬ 
tem model, MIMO radar signal model, and our radar- 
communication system spectrum sharing scenario. 

A. Coordinated Multipoint (CoMP) System 

We consider a CoMP system with a total of M base 
stations forming a set of BSs, denoted by M, and a total 
of K users forming a set of UEs, denoted by X. Each BS 
is equipped with Nbs antennas and each UE is equipped 
with A^ue antennas. The fc**' UE receives its message from 
a cluster of Mk base stations Mfc C M or the BS is 
provided with the messages of its assigned users set Xm Q X. 
Let Ufe = [uk.i---Uk,Dk]’^ C be the vector of transmit 

symbols for the UE. This complex vector represents Dk 
independent information streams intended for user k, where 
Dk < min{MkNBs, Nue). Dk is known as the degrees of 
freedom (DoF) and is defined as the number of signalling 
dimensions, where one signaling dimension is associated with 
one interference-free information stream f23i . All the BSs in 
the cluster M/c are informed about the data streams u^. This 
can be realized by utilizing the backhaul links among the BSs 
and the central switching unit. Assuming m G M^, the 
BS maps the vector via a matrix Fk,m G onto 

the signal G transmitted by the m**' BS which can 
be expressed as ll24l 

Zm — ^ ^ Fk^rn^k- (1) 

Assuming Pm to be the power constraint of the m**' BS, the 
following constraint has to be satisfied, 

E [\\im\\] = tr {E [Zmim]} (2) 

= ^ tr{Fk,mF*k^m} <Pm,rn = !,■■■ ,Mk. 
kGOCm 


The signal received by the user on the downlink is given 
by 

Yk — ^ ^ ^k,m^k,m^k ^k ( 3 ) 

where Gk^m G is the channel matrix between the 

BS and user and accounts for noise and non-radar 
interference terms. 

Assuming joint detection mode among the BSs on the uplink 
and further assuming reciprocity of the channel, the signal 
received by the m**' BS in a cluster on the uplink is given by 

bm = ^k,m^k + (4) 

k&OCm 

where G is the signal vector transmitted by the 

UE and rim accounts for noise and non-radar interference 
terms. 

B. Clustering Algorithms 

Joint transmission/reception in the CoMP system requires 
additional signalling overhead and robust backhaul channels. 
Due to this reason, only a limited number of BSs cooperate to 
form a cluster 1251. Cluster formation is an important aspect 
of the CoMP system in order to exploit benefits promised by 
the CoMP system. In general, static and dynamic clustering 
algorithms have been proposed that may or may not have 
overlapping. Each scheme has its own merits which are 
discussed as follows l2^ : 

• Static Clustering: Static clusters do not change over time 
and are designed based on the time-invariant network pa¬ 
rameters such as the geography of BSs and surroundings. 
In static BS clustering algorithms, the pre-specified BS 
clusters do not change with time and as the neighboring 
BSs interfere the most with each other on average, the 
clusters are formed by the adjacent BSs only ll25l . This 
leads to very limited performance gains as the changing 
channel conditions are not taken into account. 

• Dynamic Clustering: Dynamic clusters continuously 
adapt to changing parameters in the network such as UE 
locations and RE conditions. Dynamic clustering schemes 
exploit the effect of changing channel conditions, which 
leads to much higher performance gains without a sig¬ 
nificant overhead increase IETI, EH). So the dynamic 
clustering algorithms cluster the BSs not based on their 
geographical proximity but rather based on their level 
of interference in the absence of any cooperation 1^ . 
Thus, by capturing the effect of changing channel condi¬ 
tions, dynamic clustering algorithms exploit the inherent 
macrodiversity of multicell wireless systems [Ell. 

. Non-Overlapping vs Overlapping Clusters: Static or 
dynamic clustering algorithms can be based on schemes 
which do not support overlap, i.e. clusters are disjunct 
with respect to the cells involved. However, such a 
scheme may result in a low signal-to-interference-noise- 
ratio (SINR) for a UE at the border between clusters. This 
issue can be addressed by having spatially overlapping 
clusters. Although, overlapping clusters result in better 
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Fig. 1. Configuration of radar antenna array 


performance in terms of SINR at a UE, the price paid is 
the loss of multi-user diversity in each cluster ||26|. 

In order to keep precoder design simple and tractable we 
limit ourselves to the case where the CoMP system use static 
non-overlapping clustering algorithms where the clusters are 
formed by adjacent cells which is useful for radar since it 
will illuminate a specific geographical area. We have assumed 
that mobile network operators are performing the task of BS 
clustering and the clustering information is communicated to 
the radar during the cooperation mode. 

C. Colocated MIMO Radar 

MIMO radar is an emerging area of research and a possible 
upgrade option of legacy radar systems. Unlike the standard 
phased-array radar which transmits scaled versions of a single 
waveform, MIMO radars transmit multiple probing signals that 
can be chosen freely. This diversity in waveform enables supe¬ 
rior capabilities when compared with a standard phased-array 
radar. In this paper, we have considered a collocated MIMO 
radar with Mr antenna elements. If we denote the samples 
of baseband equivalent of M/j-dimensional transmitted radar 
signals as the signal coherence matrix can be 

written as 



1 

Pl2 

■ PlM 

i E (rz) x^j (n) = 

n—1 

... 

1 

■ P 2 M 


a; 

Pmr2 ■ 

1 


(5) 

where n is the time index, L is the total number of time 
samples and denotes the correlation coefficient between 
and signals (1 < o, c < Mr). The phases of {Poe] 
direct the beam to the angle of interest. If Poc = 0 for o ^ c, 
then Rx = Imr, i e, orthogonal waveforms. This corresponds 
to omni-directional transmission. 

The signal received from a single point target at an angle 6 
can be written as ED 

( 6 ) 


where a represents the complex path loss including the 
propagation loss and the coefficient of reflection and A{9) 
is the transmit-receive steering matrix defined as 

A(0)4at(0)a^(0). (7) 

Denoting the propagation delay between the target and the 
transmit element by t* p(0) and the delay between the target 
and the receive element by r* ;(0), the total propagation 
delay between the transmit element and the receive 
element is given by 


Tpi{9) = Tt,p{9) + Tt,i{9). (8) 


Based on these notations of propagation delays, the transmit 
steering vector, at {9) is defined as 


at{9) ^ 


,-juJcTt.l{9) p-j‘^crt,2i9) ... p-ji^cTt.MaiS) 


and the receive steering vector, ar{9) is defined as 


(9) 


ar{9) 4 


g-juiaTt,l{.6) 






T 


( 10 ) 

Using this model, the Cramer Rao Bound (CRB) for target 
direction estimation is given in equation (fTTT i lf3^ . In equation 
(fTTI i. at (9) = ^ and iir (9) = Assuming all other pa¬ 
rameters fixed, the performance of target direction estimation 
in terms of maximum likelihood or Cramer Rao bound is 
optimal for orthogonal probing signals 


D. Spectral Coexistence Scenario 

The communication system shares the spectrum with a 
monostatic ship-borne MIMO radar system as shown in 
Fig. (|2]i. The composite interference channel between radar 
transmitter and the z* cluster is denoted by ■ G 
(^MfcAfBsxMfl rpjjg superscript k implies that the cluster’s 
message is intended for UE k and 1 < z < Ct where Ct 
is the total number of clusters in the CoMP system. Channels 
are assumed to be block faded and quasi-static. The signal 
received by the z* BS cluster Mfc on the uplink, in the presence 
of radar, is given simply by 

be. = E E + + (12) 

mGCi kGOCm 

where ^xr is the interfering signal from the MIMO radar 
which we want to mitigate by designing precoders in the next 
section. 


III. Radar Precoder Design 
In this paper, according to (fT2l) . only the H matrix would 
play a role in the design of radar precoder but not the G 
or F matrices. The goal is to construct the radar precoded 
signal Scr such that its interference at all of the BSs of the 
optimally chosen cluster of the CoMP system is either zero- 
forced or minimized during interference-free mode and the 
radar keeps switching the beam across the optimal clusters 
with time. On the other hand, for cooperation mode when 
radar is communicating information to the BSs, we have to 
make sure that the bit-error-rate (BER) of the received signal 
at the BSs is low enough for effective detection. 


y(n) = aA{9)xR{n) +w{n) 
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CRB(6») 


1 

2 SNR 


(^MRa.U0)Illat{O) + at{OrRlat{0)\\ariO)f 


Mfl|an0)R^at(0)|'\ 

a*(0)Rjat(0) ) 


( 11 ) 



Fig. 2. Spectral coexistence of a MIMO radar with a coordinated multi¬ 
point (CoMP) system. An uplink-downlink model of the CoMP system while 
sharing radar bands is illustrated. 


A. Precoder Design for Interference Mitigation 

Radar interference to the i**' cluster can be mitigated 

by designing the radar precoding matrix as 

= 0, VneMfe, (13) 

= 0, VnGMfe. (14) 

The above criteria is satisfied by the following condition 

VneMfe. (15) 


Therefore, the transmitted radar signal must lie in the null 
space of ^ for all n, where n = 1, 2, • • • , Mk- We can 
rewrite equation (fTSl) as 

Pfl.^xfi e N n N . n N [Rh.R,). ( 16 ) 


Using the equality N (A) n N (B) = N (C) where C = 
[A*B*]*, the precoder must satisfy that 




(17) 


where 


H k _ 

RA — 




(18) 


To find the null space, we have to first obtain the Sin¬ 
gular Value Decomposition (SVD) of direct communication 
link iTSl which diagonalizes the channel matrix ^ G 
CMkNesxMR^ The regular form of the SVD of - is as 
follows 

= U.S.Y* (19) 


where Ui G is a unitary matrix. Si G 

(j^MfcWBsxMB ^ rectangular diagonal matrix with non¬ 
negative real numbers on the diagonal, and V* G 


is a unitary matrix. The diagonal entries of are known as 
the singular values of The MkN^s columns of Ui and 
the Mr columns of are called the left-singular vectors 
(u G and right-singular vectors (v G of 

j, respectively. 

The partitioned matrices and the outer product form of the 
SVD of j is as shown in equations (l20l i through (l22l i 
ll^ . These partitions assume that the total number of non¬ 
zero singular values, q is strictly less than M^Nss and Mr. 
From equation dTlI l. it is obvious that only the first term, i.e. 
the first q left singular vectors, Ui • • • u^, and right singular 
vectors, vi • • -v^, make any contributions to j. The rest 
of the three terms are zero. Among these three terms, the first 
term is associated with the null space corresponding to the 
zero singular values. So, for null space projection, the radar 
beam is projected onto the space formed by the second term 
of equation (ItTIi which has no contributions to j. Using 
SVD, null space of ^ is Span {Vi} where comprises 
of the columns of corresponding to zero singular values of 
As a result, the precoder PR^i must be the projection 
matrix into Span}Vi}, 

Pfl,z=V,V}. (23) 


Although this precoder will eliminate the radar interference 
to all the BSs in the i* cluster Mk, it will affect the spatial 
correlation of radar probing signals and change their coherence 
matrix. For precoded radar signals, the coherence matrix, 
Rx.i = 

Assuming that the elements of the channel matrices are 
independent and identically distributed (i.i.d.) and drawn from 
a continuous distribution, by random matrix theory, the rows 
of j are linearly independent ifJTl . The matrix Rr^ is 
therefore full rank. So, the nullity, i.e. the dimension of 
null space of ^ is greater than the difference metric, 
{Mr — M^Abs). In mathematical terms. 


null 


H k 

R,i 


= dim 


N = {Mr - MfcABs)+ . (24) 


The necessary condition for a non-trivial precoder 
(P/t.i ^ 0) to exist is that the number of radar transmit 
antennas is greater than sum of the requested DoF of all the 
base stations in a cluster. In that case, we must have. 


Mr > MkN^S 

to have a non-zero nullity for ^ and hence a non-zero 
precoder. 

In SVD, the singular values are always arranged in decreas¬ 
ing order: 

aq < aq-i < ■ • • < CTi. (25) 

These singular values are in fact measures of the trans¬ 
mitted power in the directions of their corresponding right 
singular vectors. Thus directing the radar transmitted power 
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H '' k _ 

R,i — 


Ul 


Uq+l 


UMfcAfes 



^q,(MR-q) 

0(MfcArBs-9).9 

^(MkNBS-q),{MR-q) 


Ui 


Si 



( 20 ) 


V* 


Ul 


U, 


(Tl ■ • • 0 


0 ■ • • a a 


Q J 


Ug+1 


UM^ATbs 


0 


(Mk Nss—q) R—l) 


‘'q+l 


'Mr J 


V* 


+ 


Ul 


U„ 


0 


q,{MR-q) 


'q+l 


'Mr 


Ug-j-1 • • • llMkNBS 


Ul 


U„ 


ai ■ 

• 0 


0 • 

■ 



0 


(MjjAfBs —9),9 


( 21 ) 


( 22 ) 


along vectors associated with the small singular values is 
equivalent to reducing radar interference in those directions. 
Consequently, projecting the radar transmitted signal onto a 
space spanned by the right singular vectors corresponding to 
singular values smaller than a specific threshold, crth, will min¬ 
imize radar interference to the cluster rather than eliminating 
them completely. This compromise is expected to mitigate 
the performance loss in radar target detection capability to 
some extent. The value of Uth is dependent on the power 
constraints of the communication system. With, a < ath, 
i.e. when the projection space is formed by the right singular 
vectors corresponding to the non-zero small singular values 
under a threshold in addition to the zero singular values, the 
interference would be reduced to a great extent for the current 
channel condition. This kind of projection could be called 
Small Singular Value Space Projection (SSVSP). For SSVSP, 
using SVD of j, the projection space is Span {Vs,i} where 
Vg i comprises of the columns of corresponding to the 
small non-zero singular values under ath in addition to the 
zero singular values. As a result, the precoder i is given 
by 


Pit... = Vg,,v;,. (26) 

Under switched beam condition, the radar keeps scanning 
for the best composite interference channel between itself 
and the clusters, and projects its beam towards that cluster. 
Assuming Ct represents the total number of clusters at 


a specific time in the CoMP system, with Switched Null 
Space Projection (SNSP), the best and the worst composite 
interference channels are selected in the following way iflOl 


where 


and 


where 




imax = arg max dim 
l<i<CT 




il^orst = 


ijnin = arg min dim 

l<i<CT 






(27) 

(28) 

(29) 

(30) 


Assuming full rank the nullity of ^ is defined as 


dim 


N = {Mn - Mfc7VBs)+ (31) 


where Mk is the number of BSs in the i**' cluster Mfc. 

For newly proposed Switched Small Singular Value Space 
Projection (SSSVSP), the best and the worst composite inter¬ 
ference channels are chosen as follows; 


ilhest = (H^,(32) 

where 

‘^min — arg min - Xi {||2 (33) 

l<i<CT 
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and 


where 


H„ 


= (H: 


k \. 


= arg max - xjj| 


l<i<CT 


2 • 


(34) 


(35) 


Estimation of For the simulations, we have adopted 

the channel estimation technique as described in ED. The 
estimation is carried out by having all the BSs in the i* cluster 
Mfc send training symbols to the radar and the radar utilizes 
the received signal to estimate ^ and find It is 

assumed that BSs in the cluster Mfe can cooperate in their 
choices of training symbols and transmission power. Further 
assuming channel reciprocity, the channel from BS to the 
radar transmitter is -) . The composite channel from all 

of the BSs in the cluster to the radar transmitter is given by 


H k _ 

R,i — 




(36) 


Obviously, ^ = (H^ J*. The assumption of coordination 
among BSs reduces ^ to a standard MIMO channel. So, 
by using standard MIMO channel estimation algorithm, the 
equations for estimation of composite interference channel is 
as shown below ED 


Ye = 


- - -Hi , Se 

MfciVBS 


1 < e < Lt. 


(37) 


Here, Lt is the fixed period at the beginning of each block 
of L channel uses where the estimation is performed; ye and 
We are respectively the M/j-dimensional received signal vector 
and the noise vector at time e; and p is the average SNR at each 
receiving antenna. The maximum likelihood (ML) estimation 
of the j is given by ED 


{ML) = 


I AlkN^s 


YS* (SS*) 


(38) 


where the 
W = 


Wi W2 


matrices are; Y = 
and S = 


y2 




Si S2 




SLt 


The optimal training symbols that minimizes the mean square 
error is chosen so that SS* = LtlMkNes- Cooperation among 
the BSs is required to choose the optimal training sequence. 
So, the estimation of ^ is given by 


H k 

R,i 


Hfl,. (ML) 


(39) 


B. Precoder Design for Cooperation 

In this mode, the radar communicates with the CoMP 
system. The radar informs the communication system of the 
channel matrices corresponding to base station locations that 
will be affected by its transmission and the locations that 
will have a suppressed interference due to the null space or 
small singular value space projection. The cellular network 
communicates the information about the BS clustering to 
the radar. The channel estimation is also carried out in this 
phase. But this cooperation is for brief period to take care 
of the operational security concerns of the radar. In this 


phase, radar operates in broadcast (BC) mode on downlink, 
i.e. from radar to the BSs. Simple linear precoding (LP) 
techniques such as zero forcing (ZF) and minimum mean- 
square error (MMSE) linear precoders have been considered 
for this MIMO BC mode. By inverting the channel matrix 
at the transmitter, the zero-forcing precoder can completely 
eliminate the interference among the BSs. However, a price 
of high transmission energy has to be paid especially for 
near singular matrix channels. This problem can be partially 
mitigated by using the MMSE linear precoder. This kind of 
precoder balances the transmission energy and interference 
level to achieve the minimum detection error. 

In this mode, radar signal is given by 

X]i = [xRp...XR^dR]'^ (40) 

where, dR is the number of independent information streams 
intended for the BSs and dR < min{MN bs^ Mr) and M is 
the total number of BSs in the communication system. 

Using ZE criterion, the precoding matrix is given by ll38l 

PR = ii*R(uRii*Ry' (41) 

where, H/j is the composite channel between the radar and 
all the BSs in the communication system and is given by 

[{U,R)*{n2R)*---{liMR)*]*. (42) 

The inverse in equation (l4Tli can be performed only when the 
following condition is satisfied 

Mr > dR. 

Similarly, using MMSE criterion, the precoder is as follows 

m 

= + (43) 

C. The two modes of operation and the PRI of Radar 

The radar will keep switching between these two modes of 
operation. During the first part of its Pulse Repetition Interval 
(PRI), it will operate in cooperation mode for a short period. 
But this short period is very critical for the proper operation 
of the radar for the rest of the PRI. In this brief period, it 
will collect all the information about the BS clustering of 
the CoMP and estimate the composite interference channels 
between itself and the cellular network. The radar will also 
broadcast the information it wants to communicate to the 
CoMP system. Based on the vital information it has gathered 
in the cooperation mode, it will enter next into the interference- 
mitigation mode for the rest of the PRI. Then the same cycle 
would repeat. 

IV. Spectrum Sharing Algorithms 
Based on the background theories presented in the last sec¬ 
tion about the interference-mitigation mode, we have proposed 
two new algorithms to implement Switched Small Singular 
Value Space Projection. As algorithmic implementation of 
Switched Null Space Projection is already available in litera¬ 
ture Eo), ESI, EOl, that part is omitted here. 















A. Optimal Cluster Selection Algorithm 

Our optimal cluster selection algorithm, shown in Algorithm 
O, selects the best composite interference channel based on 
the optimality criteria as mentioned in equations (l32l i and 
(|33T i. Channel State Information (CSI)s are obtained using the 
channel estimation technique as described in the last section. 
The SSVSP projection matrices of all the clusters, available at 
that time interval are then found using Algorithm (|2]l. Algo¬ 
rithm (|2]l also calculates the difference between the precoded 
and the original radar signal and returns it to Algorithm ([TJ. 
Once Algorithm ([T]i receives the difference metrics of all the 
composite interference channels, it selects the best composite 
interference channel associated with the optimal cluster and 
sends it to Algorithm (|2]i for SSVSP radar signal. 


Algorithm 1 Optimal Cluster Selection Algorithm 

loop 

for i = 1 : Ct do 
Estimate CSI of 

Send ^ to Algorithm (|2]i for small singular value 
space computation. 

Receive — x /{||2 from Algorithm (|2|. 

end for 

Find = arg min - Xi{|| 2 . 

l<i<CT 

Set H = as the best composite interference 

channel associated with the optimal cluster. 

Send H to Algorithm (|2]l to get SSVSP radar waveform. 

end loop 


B. Small Singular Value Space Projection (SSVSP) Algorithm 

In this section, we present our proposed Small Singular 
Value Space Projection algorithm. On the first cycle. Algo¬ 
rithm (|2|i gets the CSI estimates of the composite interference 
channels from Algorithm ([T]i and hnds the corresponding small 
singular value space projection matrices using the singular 
value decomposition (SVD) theorem. It also calculates the 
associated difference metrics and return them to Algorithm 
O- On the second cycle, after receiving the best composite 
interference channel matrix associated with the optimal cluster 
from Algorithm ([1]), it hnally calculates the precoded SSVSP 
radar signal by performing another SVD. 

V. Simulation Results 

In this section, the performance of MIMO radar is compared 
in terms of CRB for target direction estimation with and 
without radar precoder and as a function of number of antennas 
per BS, number of BSs per cluster, and number of radar 
antennas with either SNSP or SSVSP. We have also explored 
the effect of null space estimation error on target direction 
estimation and radar interference to clusters. The distance 
of target to radar array and the radar inter-element spacing 
are assumed to be rp = 5 km and 3A/4, respectively. The 
frequency of operation is 3.5 GHz. The signal-to-noise ratio 
between the radar and its target is denoted by SNR while 


Algorithm 2 Small Singular Value Space Projection (SSVSP) 
Algorithm 

if j received from Algorithm ([T]i then 
Perform SVD on ^ (i.e. ^ = U*S,V*). 

Find small singular value space projection matrix 

Calculate - Xfi|| 2 . 

Send - Xfl ||2 to Algorithm ([Til- 

end if 

if H received from Algorithm ([T]) then 
Perform SVD on H (i.e. H = IJSV*). 

Find small singular value space projection matrix P/j^ = 

v,v:. 

Calculate SSVSP radar signal x/j^ = P/j^x/j. 

end if 



Fig. 3. CRB on direction estimation of the target as a function of antennas 
employed by the BS when using NSP and SSVSP in a spectrum sharing 
scenario between MIMO radar and cellular system. A cluster size of = 3 
is used while having Mr = 100 radar antenna elements. 

the same parameter between the radar and the cluster of BSs 
of CoMP system is represented by p. The target direction is 
assumed to be at 0 = 0°. 

A. Performance Analysis of Interference Mitigating Precoder 

Fig. (|3l shows the radar CRB performance for orthogonal 
radar signals as well as precoded radar signals for NSP and 
SSVSP with estimated Channel State Information (CSI). The 
figure reflects the fact that SSVSP performs better than NSP 
from the perspective of radar target detection capability as 
predicted. This plot also shows that as Wbs increases, the target 
localization performance of radar degrades. Because with 
the increase in Wbs, {Mu — MkN^s) decreases for constant 
Mfi and M^- So, null space of ^ shrinks with consequent 
impact on the precoder which in turn degrades CRB(0). 

In Fig. dUl, we have investigated the effect of number of 
radar antennas {Mji) on the CRB performance of radar. As 
obvious from equation (ITTIl that CRB explicitly depends on 
Mu and improves when M/j increases. For precoded radar 
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Fig. 4. Effect of number of radar antennas (Mr) on the CRB(0) performance 
of radar. A cluster size of = 3 is used and N^s = 8. 



Fig. 6. Radai' Interference on one BS with number of radar antenna elements, 
M^. Abs = 8 antenna elements per BS are used while having cluster size, 
Affc = 3. 




Fig. 5. Variation of radar interference on BS clusters with cluster size (M/.). 
TVbs = 6 BS antenna elements are used while having Mr = 100 radar 
antenna elements. 


signals, increase in Mr leads to increase in the nullity of ^ 
which impacts the choice of and results in improved CRB 
performance. Results in this plot indicate that increasing the 
number of radar antennas can compensate for the performance 
degradation in target direction estimation due to correlation in 
the precoded radar signals. Consequently, for a given target 
RMSE, number of radar antennas must increase when radar 
employs a precoder to zero-force or minimize its interference 
at the clusters. 

In Fig. dU, the effect of interference at the clusters is mea¬ 
sured using the metric, ^nR ■ Interference 

becomes less significant as the duration of training phase, Lt 
increases while it will be completely eliminated with perfect 
knowledge of So, estimated CSI approaches perfect CSI 
as Lt increases. This plot also reflects the fact that from the 
perspective of interference at the clusters, NSP outperforms 
SSVSP to a great extent. Fig. (|6]l represents the same facts 
for a specific BS in a cluster. Fig. (|7]i plots radar target 


Fig. 7. CRB on direction estimation of the target as a function of number 
of antenna elements employed by the radar (Mr) with switched NSP and 
SSVSP while in a spectrum sharing scenario between MIMO radar and 
cellular system. A cluster size of Mk = 3 is used while having estimated 
CSI. Radar beam is swept across 4 clusters (Cr = 4) and the nullity of 
the clusters is varied by varying the number of antenna elements per BS, 
[Afss = (6, 5, 4,3)]. 


estimation capability with switched NSP and SSVSP. It is 
obvious that both of the beam sweeping methods improve 
radar performance when the beam is projected on to the 
optimal cluster associated either with the maximum null space 
or minimum deviation of the precoded radar signal from the 
original one. For the worst-case scenarios, SSVSP outperforms 
SNSP 

B. Performance Analysis of Information Exchange Precoder 

Fig. ([8]l through Fig. (fTOl l show the spectrum sharing sce¬ 
nario during cooperation mode. Fig. ([8]) and Fig. (|9]l show 
that target location performance of both ZF and MMSF 
precoding is same at high p for low values of TVbs while they 
start deviating from each other at high TVbs. ZF precoding 








































10 





—B— ZF precoder (estimated CSI, Mj^=30) 

0 MMSE precoder (estimated CSI, Mp|=30) 
—«— ZF precoder (estimated CSI, Mj^=100) 

— 1 — MMSE precoder (estimated CSI, Mp|=100) 




0 5 10 15 

SNR between radar and comm system, p(dB) 


Fig. 8. CRB on direction estimation of the target as a function of the number 
of antenna elements employed by the BS while in a cooperative spectrum 
sharing scenario between MIMO radar and cellular system. A total of 3 BSs 
(M = 3) is used while having estimated CSI and Mji = 26 radar antenna 
elements. 



Fig. 9. CRB on direction estimation of the target as a function of the number 
of antenna elements employed by the BS while in a cooperative spectrum 
sharing scenario between MIMO radar and cellular system. A total of 3 BSs 
(M = 3) is used while having estimated CSI and = 100 radar antenna 
elements. 

outperforms MMSE precoding at low p. They also reflect the 
fact that this discrepancy in performance can be overcome by 
increasing the number of radar antenna elements. Fig. (Unii 
represents the performance of the cluster in terms of BER. 
As predicted before, ZF outperforms MMSE precoding in 
this case too, but at the expense of high transmission energy. 
The increase in the number of radar antenna is the obvious 
solution to close the gap between the performance of these 
two precoding methods. 

VI. Conclusion 

In this paper, the authors have designed precoders of a 
radar in a scenario where a MIMO radar co-exists spectrally 
with a CoMP commercial communication system. We have 
investigated the design of precoder for two modes of radar 


Fig. 10. BER performance during cooperation mode using QPSK waveform. 
A total of 4 BSs {M = 4) is used. 

operation: interference-mitigation mode, when the radar at¬ 
tempts to avoid interference with the communication system 
and cooperation mode, when it exchanges information with 
the CoMP system. For the interference-mitigation mode, the 
radar steers its beam towards the optimal clusters of base 
stations and keeps switching the beam from one optimal 
cluster to another. In case of traditional SNSP, the optimality is 
decided upon the maximum nullity of the cluster. A new space 
projection and switching method has been proposed where the 
beam is steered to the small singular value space of the optimal 
cluster whose optimality is chosen based on the minimum 
difference between the precoded and original radar signal. 
For cooperation mode, the precoder is designed based on the 
conventional zero- forcing and MMSE criteria to reduce BER 
at the BSs for effective detection. Channel estimation error 
has been taken into account for both modes of operation. Al¬ 
though the introduction of precoding tries to mitigate the radar 
interference at the clusters during interference-free mode and 
minimize interference among BSs during cooperation mode, 
in doing so it causes the radar probing signals loose their 
orthogonality with consequent performance loss in radar’s 
target parameter estimation capability. SSSVSP mitigates that 
effect to some extent but increase in the number of radar 
antenna elements has stronger impact. 
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